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Abstract 

We construct and solve the boundary Yang-Baxter equation in the RSOS/SOS 
representation. We find two classes of trigonometric solutions; diagonal and non- 
diagonal. As a lattice model, these two classes of solutions correspond to RSOS/SOS 
models with fixed and free boundary spins respectively. Applied to (l+l)-dimenional 
quantum field theory, these solutions give the boundary scattering amplitudes of 
the particles. For the diagonal solution, we propose an algebraic Bethe ansatz 
method to diagonalize the SOS-type transfer matrix with boundary and obtain 
the Bethe ansatz equations. 
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1 Introduction 



In the study of the two-dimensional integrable models of quantum field theories and 
statistical models, the Yang-Baxter equation (YBE) plays essential roles in establishing 
the integrability and solving the models. In the field theories, the YBE provides a 
consistency condition for the two-body scattering amplitudes (5-matrices) in the multi- 
particle scattering processes since the scattering is factorizable. With unitarity and 
crossing symmetry, the YBE can determine the S'-matrix completely, although not 
uniquely due to the CDD factor which is any function of the rapidity satisfying the 
unitarity and crossing symmetry conditions. This CDD factor is usually neglected under 
the the minimality assumption. In the statistical models, if the Boltzmann weights 
satisfy the YBE, row-to-row transfer matrices with different values of the spectral 
parameter commute each other so that the models are integrable. 

Recently there has been a lot of efforts in extending these approaches to models with 
boundaries. The main motivation is that these models can be applied to 3D spherically 
symmetric physical systems where s-wave element becomes dominant. One-channel 
Kondo problem, monopole-catalyzed proton decay are frequently cited examples. Also 
one can generalize the conventional periodic boundary condition of the statistical mod- 
els to other types like the fixed and free conditions. 

The existence of the boundary adds new quantities like boundary scattering ampli- 
tudes and boundary Boltzmann weights, and one needs to extend the YBE to include 
these objects. The boundary Yang-Baxter equation (BYBE) (also known as the re- 
flection equation) |l]] plays the role of the YBE for the integrable statistical models 
P, |3| and quantum field theories [Q] in the presence of a boundary; it is the necessary 
condition for the integrability of these models. The equation takes the form 

Ri(u)Si 2 (u +u)R 2 (u)Si 2 (u - u) = Su(u - u)S 2 (u )S\ 2 (u +u)Ri(u) (1.1) 

where R\( 2 ) is the boundary scattering (or reflection) matrix in the auxiliary space 1(2) 
and S12 is the solution to the YBE. In general, R(u) need not be a C-number matrix, so 
the equation may be taken as the defining relation for the associative algebra generated 
by the symbols R(u) ||. This algebra possesses a very rich structure and has been found 
to be connected with braid groups 0], lattice current algebra M, twisted Yangian Q 
and so on. Taking the quantum space of R(u) to be trivial, the BYBE is a quadratic 
matrix equation which allows the matrix R(u) to be solved for given S(u). To date, 
several solutions of the BYBE have appeared in the form of vertex representation, far 
less is known however for the solution in the solid on solid (SOS) or restricted solid on 
solid (RSOS) representation [||, pLQ] . Also less clear is the vertex-SOS correspondence 
associated with this algebra. In particular there is no reason to expect the vertex-SOS 
transformation for the YBE continue to hold for this algebra. Therefore finding solution 
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in the RSOS form may help to clarify the issue of the vertex-SOS correspondence. 
Moreover, the RSOS solution will reveal the special mathematical structure associated 
with this algebra when the deformation parameter is a root of unity. 

From a physical point of view, the solutions have applications in statistical me- 
chanics and field theory. In the context of statistical mechanics, the solutions give rise 
to integrable SOS/RSOS models with boundaries where the C-number solution of the 
BYBE provides the Boltzmann weights of the statistical models at the boundary. The 
first nontrivial case gives the tri-critical Ising model. The study of integrable statistical 
model with boundary will shed light on the issue of the dependence of the Casimir 
energy on the boundary and surface properties |11, 27]. From the field theory point of 
view, the solutions are relevant to the study of the restricted sine-Gordon model [16] 
and the perturbed (coset |]l^, [l3|]) conformal field theory (CFT) [14] with boundary. 
In this case, solutions to the BYBE are scattering matrices of the particles with the 
boundary. 

In this paper we construct the BYBE in the RSOS/SOS representation. The equa- 
tion is studied for the diagonal and non-diagonal cases, and the most general trigono- 
metric solutions are found up to an overall factor. This factor is then fixed using the 
boundary crossing and unitarity conditions § (up to the usual CDD ambiguity). We 
then construct an integrable RSOS/SOS model using the diagonal solution and propose 
an algebraic Bethe ansatz method to diagonalize the transfer matrix of the SOS model. 



2 Solutions to the boundary Yang-Baxter equation 
2.1 Generalities 

In this section we solve the BYBE for the RSOS(p) ; p = 3, 4 . . . scattering theory. The 
RSOS(p) scattering theory is based on a (p — 1) - fold degenerate vacuum structure, 
vacua can be associated with nodes of the A p -\ Dynkin diagram. The quasi particles in 
the scattering theory are kinks that interpolate neighboring vacua, they can be denoted 
by non-commutative symbols K ab {u) where \a — b\ = 1 with a, b = 1, . . . ,p — 1 and u 
is related to the the kink rapidity 6 by u = —i9/p, so that the physical strip is given 
by < Re u < it/ p. In the rest of the paper, we will refer to a, b as heights or spins. 
Formally, scattering between two kinks can be represented by the following equation 
(see Fig.©) 

K da {u)K ab {u) = Y,Sf c {u - u')K dc ( U ')K cb (u) (2.1.1) 
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where the S'-matrix is given by 



St(u)=U(u)(^j U ' 2l W£{u) (2.1.2) 



and 

1/2 



Wg(u) = Lmu5 bd (j*) + sin( 7 - u)8 ac I (2.1.3) 



satisfies the YBE in the RSOS representation. 

b 




d 

Figure 1: The bulk RSOS scattering matrix S°£(u — u) 

Here [a] denotes the usual (/-number given by 

r , sin(<27) 7r 

[a] = — 7 = - 

sin 7 p 

and the overall factor U{u) is a product of Gamma functions satisfying the relations 

U(u)U(— u) sin(7 — u) sin(7 + u) = 1 

W(7 -u) = U(u) , 

and is given by 



vr ttV^ F,(0)F,(7) ' 



r ( ?tL _ tA r ( i + ^ _ A 

- / (2i+l) 7 tt \ r / t (2f-l) 7 »\ • 
\ 7T n J \ it it J 
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This factor, together with the overall ^-number factor, ensures that the S'-matrix sat- 
isfies both crossing and unitarity constraints: 

S&(u) = S b a c d {l-u) (2.1.5) 
Y^Sf,{u)si h {-u) = 5 ac . (2.1.6) 

c 

Consider now the above scattering theory in the presence of a boundary denoted 
formally by B a , then the scattering between the kink and the boundary is described 
by the equation 

K ab (u)B a = R b ac (u)K bc (-u)B c (2.1.7) 

c 

which can be given a graphical representation shown in Fig.(j2|). Notice that in this 
representation, the boundary naturally carries an RSOS spin. 

b 



a 




c 



Figure 2: The boundary RSOS scattering matrix R b ac (u) 

The function R b ac is called the boundary scattering matrix and satisfies the BYBE, 
which in the RSOS representation takes the form 

]T R a bb , (u)S a / a , (u + u)R a b ; b „ {u)Sy a „ (u -u) = 

a ,o 

£ S£ (u - u)R a bb , {u)S«% (u + u)Rf b „ («) . (2.1.8) 

' if 

a ,o 

This equation is illustrated graphically in Fig.(||). 
In general, the function R bc (u) can be written as 

RUu) = n{u) iiU [h^XUv) + he {5 b , a+l U a {u) + h,a-lDa(u)}] (2.1.9) 

\[a\[a\J 

where lZ(u) has to be determined from the boundary crossing and unitarity constraints, 
while X bc and U a , D a have to be determined from the BYBE. An overall q- number factor 
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b 



b' 



b" b 



b' b' 



Figure 3: The boundary Yang-Baxter equation 

has also been multiplied to the above to cancel that from the bulk S'-matrix in order 
to simplify the BYBE. If X% c does not vanish, the boundary i?-matrix describes non- 
diagonal scattering process, otherwise the scattering is called diagonal. Note that due 
to the restriction that vacuum assumes value 1, ... ,p — 1, X^ c , X£~ , D\, U p -\ are not 
defined. The case p = 3 has only diagonal scattering, so X£ c does not exist. 

2.2 Non-diagonal scattering 

We consider the scattering where the off-diagonal component X% c is non-vanishing. To 
start, consider the case b ^ c ^ b in eqn. fl2,l!8| ) where the BYBE gives 

X2-x, a+ i(u')X:+l a+3 (u) = X a a ^ a+1 {u)X a a Xl a ^{u) ; 2 < a < p - 4 (2.2.1) 

which implies that can be written as 

X a a±w {u) = h ± {u)X a ± (2.2.2) 

where h±(u) depends only on u and X± only on a. 

On the other hand, the case c = b = b" , a = a" gives 

XZ-l,a+l(u')XZ + l, a -l(u) = *£-l,a + l(«)*a+l,a-l(«') ^ < a < p - 2 (2.2.3) 
which implies that 

h+{u)h-{u) = h + {u)h^{u) , (2.2.4) 
from which we conclude that 

h + (u) = (const. )h-(u) ■ (2.2.5) 

Absorbing the constant in the above equation into X° or X!t, we can make h + equal 
to h- so that we can absorb the h±{u) into the overall TZ(u) factor and treat Xg c as u 
independent from now on. 
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With this simplification, eqn.(2.1.8) can be broken down into the following inde- 
pendent equations in addition to the above two equations: 



U a D a+2 f+ 1 + /. 



[a+1] 



+ D a+2 D a+2 f_ [1 + f. 



+ -^o+l,a+3^a+3,a+l /- - U aD a+2 f+ I 1 + /- 



[a + 2] 
[a + 1] 
[a + 2] 
[a + 1] 



+ U a U a f- (l + /+ — J + ^0-1,0+1^0+1,0-1 / 



for 1 < a < p — 3, 



D a+1 f- 1 + f- 



[a+1] 



+ [/._!/+ ! + /_ 



[a-1] 



C/a-l/+ " C/a+l/- 



(2.2.6) 



(2.2.7) 



UJ-(l + f. 



[a + 1] 
D a+2 f+ - DJ. 



+ D a+2 f+ l + /_ 



[a + 2] 
[a + 1] 



for 2 < a < p — 3, and 



C 2 /+/-^ [ _r I pr-^ + ^o(i + /- 

= f 1 + / 4 



[a-1] 



[a-1] 
C/o(l + /- ' 



1 + /-+ 



[a-1] 



[a-1] 



= u a (i + u 



[a + 1] 



1 + /4 



[a + 1] 



[a + 1] 



D a 1 + /_ 



[a + 1] 



(2.2.J 



(2.2.9) 



(2.2.10) 



for 2 < a < p — 2. The last four equations are derived based on the assumption that the 
off-diagonal weight X% c is nonvanishing. In the above equations, we used a compact 
notation where U a = U a (u), U a = U a (u ) (similarly for D a ) and 

/+ = sin(ii ± u)j sin(7 — u + u) . 

In addition, it should also be mentioned that the last term in the rhs (lhs) of 
eqn.( [2.2."6| ) is present only when a ^ l(p — 3) and the first terms of eqns.( 2T2~lj| ) and 
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( 2.2.10 ) are allowed only for a / 2 and a ^ p — 2, respectively. Let us call these 
terms that are not supposed to be there the "unwanted" terms. The way to solve these 
equations is to construct some recursion relations for the functions U a , D a , X®±i aZfl 
and solve the recursion relations subject to the conditions that the above mentioned 
"unwanted" terms are zero. Explicitly, one has 



Y 1 Y 1 

^0,2-^2,0 



yP-l yP _1 

^-p-2,p^p,p-2 
= D n 







(2.2.11) 
(2.2.12) 



m 



A few comments are in order here. Notice that the coefficients of Uq and D p 
eqns. (|2.2.9D , ( [2.2.10| ) are [0] and \p] respectively, which vanish by construction. So 
in principle, one needs not impose the above condition, eqn.( 2?27L^ ), for the recursion 
relations of U a ,D a . However, we shall see later that a particular solution of U a ,D a is 
given by 

U a = -D a = (const.) (2.2.13) 
which cancels the vanishing coefficients at a = 0,p, and renders the unwanted terms 



nonvanishing. Therefore, we have to impose eqn.(2.2.12) on the above solution. 

Notice also that most of the above equations do not apply to the case p = 4, so we 
shall deal with this case separately. 

From the above it is clear that eqns.( 2T2^ )-( |2.2!Tl| ) can be divided into two parts; 
eqn.( |2.2.8D -eqn.( |2.2.1C| ) determine U a and D a while eqns. (|2.2.6|) ,( p.2.11| ) determine Xj} c . 
Indeed, comparing eqn. fl2~2~7D with eqn.( l2~2~9| ) and similarly eqn. (fE2l| with eqn.( |2.2.iq ), 
we deduce that 



(D a -D a ) 1 + / 



[a-1] 



(U a _ 2 - C/ a _ 2 )/+ 

+ (U' a ~ U a ) 



[a-1] 



(U a -U a )[l + f. 



[a + 1] 

Substituting one into another, we get 

(U' a+2 - U a+2 ) - (U' a - Ug) 
{U' a - U a ) - (U' a _ 2 - U a -2) 

(D' a+2 -D a+2 )-(D' a -D a ) 



^a+2)/4 



a+2 „a^,^ [a + l] 
+ (D' a -D a ). 



sin {^(a + 1)7 + u + u 
sin ((a — 1)7 + u + u 
sin ( (a + 1)7 — u — u 



(D' a - D a ) - (D' a _ 2 - D a - 2 ) sin ((a - 1)7 - u 



(2.2.14) 
(2.2.15) 

(2.2.16) 
(2.2.17) 
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and writing the rhs respectively as 



cos hu + (a + 1)7) - cos (2u + (a + 1)7) 
cos ( 2u + (a — 1)7 ) — cos (2u + (a — 1)7) 



cos (2u — (a + 1)7^) — cos (2u — (a + 1)7) 
cos ( 2u — (a — 1)7 ) — cos (2u — (a — 1)7) 



it is clear that 

U a +2(u) - U a (u) = cos (2u + (a + 1)7) + fi' a 
D a+2 {u) - D a (u) = - cos (2u - (a + 1)7) + (f> a 

where a , (f> a are unknown functions of a only. Iterating the above, one finds 

U a {u) oc sin(2u + 07) + a(u) + (3 a (2.2.18) 
D a (u) oc sin(2u - 07) + j(u) + (f> a (2.2.19) 

where a(u),^(u) are unknown functions of u, and f3 a ,4>a of a. Furthermore, from 
eqns. (p^2~^ )-( |2.2!To|) , one can establish the following symmetry 

U a {u) = -D a {-u) , (2.2.20) 

which reduces the number of unknown functions to two, namely, a(u) and j3 a - To 
determine them, we have to substitute the above expressions for U a ,D a back into 
eqns.( 2^2~7| )- (|2.2.10D . Notice, however, that these equations are linear in U a ,D a , so it 
suffices to consider a(u) and (5 a separately. Doing this amounts to finding special solu- 
tions to eqns. ( |2.2.7D - (|2.2.10D where U a ,D a have only u or a dependence. The solutions 
are given by 

a(u) = (const.) . ] , (3 a = (const.) — ^ — - (2.2.21) 

sm(2u) sm(<27) 

respectively. Imposing eqn.( |2.2.T2D on a(u),(3 a , one finds 

P a = . (2.2.22) 

So far the above U a , D a are obtained based on the assumption that the numerators and 
denominators of the lhs of eqns. ( ^72 . 1 6| ) , ( ^72 . 1 7| ) are nonvanishing. In fact, there are 
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additional solution to these recursion relations where the above mentioned numerators 
and denominators vanish; 

U a+2 {u)-U a {u)-U a+2 {u) + U a {u) = 
D a+2 {u)-D a {u)-D a+2 {u)+D a {u) = 0. 

The additional solution is given by 

. . , , sin(a7 + 2u) . . , . sin (try — 2u) 

U a {u) = const. . ,\ { . \ , D a {u)= const. . ( \ { . . \ ■ 2.2.23 
sm(2u) sm(a7j sm(2n) sm(a7) 



However, on imposing the condition eqn.( 2.2.12 ), one finds that they have to be zero. 
In summary, the general non-diagonal solution to the four linear equations is given by 

U a (u) = sin(2n + aj) + (const.) . L , 

sm{Zu) (2.2.24) 
D a +i(u) = sin (2u - (a + 1)7) + (const.) sin ( 2n ) 

where 1 < a < p — 2. 

Having found U a ,D a , the function X^ c can be easily obtained from eqn.( 2.2.6| ), 
which can be further simplified with the symmetry properties given in eqn. ( 2.2.20 ) and 
taking u to be — u since Xg c does not depend on the rapidity. This gives 

Ktla+sKtla+l ~ K-l,a+lK + l,a-l = D a+2 (u)U a+2 ( U ) - D a {u)U a {u) . 

Substituting U a , D a into the rhs and iterating the equations, we get 

K^a+iK+i,a-i = sm 2 7 - sin 2 (a 7 ) (2.2.25) 

where the use of eqn. ( |2.2.11 ) forces the first term to be sin 2 7 and the coefficient of sin ^ 2u ^ 
to vanish. Since this equation determines only the product, a+1 and X^ +1 a _ 1 are 

determined up to a gauge factor. 

These solutions have the property that 

Up- a (u) = —D a (u) , X p _®_ lp _ a+1 Xp_2 +lp _ a _i = X^_ l a+1 X^ +l a _ 1 . (2.2.26) 

Next, we consider the BYBE for p = 4. The functions U 2 ,D 2 and X^^X^ satisfy 
the following equations 

xUu)Xl{u) = XUu)Xl(u) (2.2.27) 
U 2 (l + V2/_) + D' 2 (l + y/2f+) (l + y/2f. 



= U 2 + D 2 [l + V2f + ) (2.2.28) 
D 2 (l + V2f_) +U' 2 (l + V2f+) (l + V2/- 



D' 2 + U 2 (l + V2f+) ■ (2.2.29) 
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The functions Ui,D$ are diagonal scattering components and do not couple to the 
above equations and we shall defer to next section for their computation. Here we have 
used the compact notation introduced earlier for U a , D a , and written out explicitly the 
rapidity dependence of Xg c . As before, the last two equations have been derived based 
on the assumption that Xf s ,X 31 are nonvanishing. From eqn. fl2,2l2^ ), we deduce that 

Xf 3 (u) = (const. )Xl l {u) . (2.2.30) 

We are free to take X 31 as unity and the above implies that Xf s is just a gauge factor, 
which we call g. 

The rest of the equations can be turned into ordinary differential equations in the 
limit u — > u, giving 

(u 2 (u)+D 2 (u)^tan(2u) + 2(U 2 (u)+D 2 (u)) = (2.2.31) 

(u 2 (u)-D 2 (u))cot(2u)-2(U 2 (u)-D 2 (u)) = 0, (2.2.32) 

which can be integrated to give 

U 2 (u) = B / sm(2u) + C cos(2u) (2.2.33) 
D 2 (u) = B / sm(2u) - C cos(2u) , (2.2.34) 

with B, C as free parameters. 

This completes the determination of the non-diagonal solutions of the BYBE. 

2.3 Diagonal scattering 

For the diagonal scattering, we take 

RUu) = {[b]/[a))- u ^ll(u)5 bc [5 bia+1 U a (u) + 5 b ^ x D a {u)\ , (2.3.1) 
and the BYBE is equivalent to a single equation 

U a -i (u )D a+ i (u) sin(u + u) sin(a7 — u + u) 

+ D a+ i(u )D a+ i(u) sin(w - u) sin(a7 + u + u) = 

U a -i{u)D a+ i(u ) sin(n + u) sin(a7 + u — u) 

+ U a -i(u )U a -i{u) sm(u — u) sin(a7 — u — u) , (2.3.2) 

which holds only for 2 < a < p—2. So the functions U p - 2 and D 2 can not be determined 
from the BYBE. The above can be recast into 

[cos(a7 + 2u)Z(u) — cos(a7 — 2u)][Z(u ) — 1] = 

[cos(a7 + 2u)Z(u) - cos(a7 - 2u)][Z{u) - 1] (2.3.3) 
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where 

Z a (u) = D a+ x{u)/U a -.\{u) . 
One can easily find that the general solution is 

Dg+i{u) = sin(g a + u) sin(g a + aj - u) ^ g 

U a -i(u) sin(£ a - u) sin(£ a + 07 + u) 

where £ a is a free parameter. 

Thus for the diagonal solution, there are p — 3 parameters £ a . This solution gives 
p — 1 distinct diagonal scattering theories of kinks, each with a specific boundary B a . 
There is one free parameter £ a for each theory, except for the cases of a = l,p—l where 
there is no free parameter. 

This solution includes a particular case of p = 4 which has been omitted earlier. 

Further relations from boundary unitarity and crossing symmetry will be required 
to disentangle U a -\ and D a+ i, and determine JJ% and -Dp-2, see later. 



2.4 Boundary unitarity and crossing symmetry 

The boundary unitarity and crossing symmetry conditions of the scattering matrix 
R^ c {u) determine to some extent the overall factor R(u). These conditions can be 
written as 

J2RUu)R a cd (-u) = 5 bd (2.4.1) 

c 

J2S%(2u)R d bc ( 7 /2 + u) = RUl/2-u). (2.4.2) 
d 

Consider the non-diagonal scattering (p > 4) first. Substituting the expression for 
Rfc into the unitarity condition, we get the following 

K(u)n{-u) [W^ - U a (u)D a (u)] = 1 

where use has been made of the symmetry eqn. (I2.2.2O0 . Applying the results eqns.( |2~2.24D , 
(2.2.25) to the above leads to 

n{u)n{-u) (sin 2 7-4 sin 2 u + 4 sin 4 u) = 1 . (2.4.3) 

While for crossing symmetry condition we get 

U(2u)TZ('j/2 + u) sin( 7 - 2u) = K(j/2 - u) (2.4.4) 
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using the relations 



[o+l] " w > + l] 
= f{2u)(U a (u)-U a { 1 -u)) 



+ -D a +2(7 - u) 



[a + 2} 



[a + 1] a ^ yi > + l] 
= f(2u) (D a+2 (u) - D a+2 ( 7 - u)) 



(2.4.5) 



(2.4.6) 



which are obtained from eqn.( 2.2!8|) in the limit u — > 7 — u. Here f{u) = sinu/ sin(7 — 
u). 

The factor TZ(u) can be determined from eqns. (|2.4.3|) , (I2.4.4D up to the usual CDD 
ambiguity by separating lZ(u) = 1Zq{u)TZ\(u) where TZq satisfies 



K (u)K (-u) = 1 

W(2n)^ (7/2 + n)sin(7-2n) = ^0(7/2 -it) 



whose minimal solution reads 



Ra(u) 



Fq(u) 
F (-u) 



where Fi(u) has been given in eqn,( |2.ll4 ). While 1Z\ satisfies 

TZi (u)lZi(-u) (sin 2 7-4 sin 2 u + 4 sin 4 u) = 1 
Ki(u) =^i( 7 -u) 



with minimum solution 



(2.4.7) 
(2.4.8) 

(2.4.9) 
(2.4.10) 



where 

<r(x, u) 

n(x,«) 



no*, i - u) n(-*, i-u) u(x, -i+u) u(-x, -1 + u) 



lFM)if(-x, i) 



n 



r(l + (2Z + l)^ + f -f 

or(i + (2J + §)i + f-* 



(2.4.11) 



For p = 4, the crossing symmetry condition is the same as eqn.(2A4|), but unitarity 
now requires that 



K(u)K(-u) {g + C 2 cos 2 (2u) - B 2 / sin^ 



(2.4.12) 
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These equations can again be solved by separation (see [17] for details). It should be 
remarked that there are two parameters in this case while there is only one in the higher 
p cases. 

Next, we consider the diagonal case. For convenience, we set the undefined terms 
Uo,Dp to be zero. Unitarity and crossing symmetries relations give respectively 



n(u)n{-u)u a {u)u a {-u) = 1 

K(u)K(-u)D a+1 (u)D a+1 (-u) = 1 



(2.4.13) 



and 



U(d — 2u)1Z{^ — u) [f7 (7 — u) sui7sin(a7 + 2u) 
+D a+2 {l - it) sin( 7 - 2u) sin ((a + 2)7)] = K{u)U a {u) sin((a + 1)7) 

U(^i — 2u)lZ{^ — u) [D a+ i (7 — u) sin 7 sin ((a + 1)7 — 2u) 
+U a -i{l ~ u) sin(7 - 2u) sin ((a - 1)7)] = TZ(u)D a+1 (u) sin(a7) 

for 1 < a < p — 2. 

These equations can be solved separately; we set 

K(u)K(-u) = 1 

U(2u)TZ(-f/2 + u) sin(7 - 2u) = K(j/2 - u) , 



(2.4.14) 



(2.4.15) 



so that TZ(u) has exactly the same solution as that of TZo(u) considered earlier. While 
U a , D a satisfy 

U a (u)U a (-u) = 1 



D a+ i(u)D a+1 (-u) = 1 



and 



U a (j - u) sin7sin(a7 + 2u) + D a+2 (j - u) sin(j - 2u) sin((a + 2)7) 

= U a (u) sin(2n) sin((a + 1)7) , 
Ai+i(7 _ u ) sin7sin((a + 1)7 - 2u) + J7 a -i(7 - u) sin(j - 2u) sin((a - 1)7) 

= D a+ i(u) sin(2u) sin(a7) . 

(2.4.17) 

for 1 < a < p — 2. Substituting eqn.( p.3^ ) into the above we get a relation between 
U a (u) (£>«(«)) and U a (j - u) {D a (j - «)); 

Ug-i(u) sin(2(7 - u)) sin(£ a - u) sin(^ a + 07 + u) 

U a -i{^-u) sm{2u) sin(£ a - (7 - u)) sin(£ a + 07 + (7 - u)) 

(2.4.18) 

Da + i{u) sin(2(7 - u)) sin(g a + u) sin(g a + aj - u) 

D a +i{l ~ u) sin(2u) sin(£ a + (7 - u)) sin(£ a + 07 - (7 - u)) 
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These relations together with eqn.( |2.4.1~6 ) give the minimal solutions of U a (u) and 
D a (u) 

U a -i{u) = ^ sin 2(7 - ii)sin(£ a - u) sin(£ a + £17 + u)a(j,u) 

x - 7, «)o-(5 - &» - - 07, «) , (2.4.19) 

D a +l(w) = ^sin2(7-ti)sin(^ a + u)sin(^ a + a7-u)o-(7,u) 

x - 7, uM| - fa, «M| - - 07, «) • (2.4.20) 

for 2 < a < p — 2 and 

£W«) = ^sui2(7-«)(t(7,«M|-7,«), (2.4.21) 

D 2 (u) = isin2( 7 -n)a(7,n)a(| - 7 ,n) . (2.4.22) 

To summarize, there are two classes of solutions to the BYBE; diagonal and non- 
diagonal. Unlike the solution in vertex representation, the former is not a special 
limit of the latter. In fact, the diagonal solution carries p — 3 free parameters while 
nondiagonal has non. 

2.5 Comments on the SOS model 

We have considered from the start heights take values from 1 to p— 1, which is necessary 
for the bulk scattering weights to be finite as the parameter tt/j = p is a positive 
integer. When ir/"f is not a rational number, there is no bounds on the heights and 
the corresponding representation is known as solid-on-solid (SOS). The removal of the 
heights' restriction certainly affects the solution of the BYBE. 

For the diagonal solution, it is clear that essentially there is no difference between 
the SOS and RSOS solution, where the SOS solution is given by eqns.( 2~4.19| ),( 2~4.20D 



for any integer a and the overall factor lZ{u) is the same as before. Hence there is a 
free parameter for each height. 

While for the non-diagonal solution, the conditions eqns. (12.2.11]) , (|2.2.12| ) of the 
recursion relations do not apply at all, thus /3 a 7^ and the corresponding diagonal 
weights are given by 



tt 1 \ ■ to 1 \ 1 1 1 sm(a7 + 2u) 1 

U a (u) = Ci sm(2u + 07) + 0> — — j — \ + C3— — \ . — 7 — -x + Ca . jn \ 

aK ' 1 v " z sm(a7) d sm(2n) sm(a7) 4 sm(2it) 

D a (u) = ci sin(2ii — aj) — Co— — \ — r + C3 ?^{ a 7 . n + ca . L N 

ay ' y " z sm(a7) d sm(2n) sm(a7) 4 sm(2u) 



(2.5.1) 
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where Cj ; i = 1, . . . , 4 are free parameters. The off-diagonal weights satisfy 
X a-l,a+l X a+l,a-l = 4 + 2c i c 4 cos(a 7 ) - cf sin 2 (a 7 ) - c| I 



sin (07) 



cos(a7) 9 1 , v 

- 2c 2C 3-^- - ci-^-— (2.5.2) 
sm (07) sin (07) 

and are determined up an additive constant cq. Of these five free parameters one is an 
overall factor, so there are four free parameters in this case. The overall factor TZq is 
given as before, but 1Zi(u) now contains all the information of the boundary conditions 
and has to satisfy 

ll^K^-u) (c\-c\ sin 2 (2n) + 2 Cl c 2 cos(2u) - 4 + ^ 3 c A cos(2u) + c 2 A = 1 

V sin (2u) / 

K 1 {u)=K 1 (~f-u) 



(2.5.3) 



whose minimum solution is 



(c 3 + C 4 )cr(0,tt) 



where are related to c,- via 



c 3 + c 4 x 2 
4ci 



Yl cos 2 tfl = 
1=1 

2 q 2 Q 2 q c + C 3 C 4 + 2ciC 2 

> COS Vi COS U,- COS Vk = o , 

-t<J<fc=l 1 

V cos 2 cos 2 i?j = C ° + 5C1C 2 2 + 4C ' , 
J<J=1 1 

£cos 2 ^ = 

i=l c l 

and <r(x, u) has been given before, while 

c f?) _ft r (^ + 1 )^ + f) r (^ + 2 )^ + f) r ( 1 -^ + 1 )^-f) r ( 1 -^ + 2 )^-f) 
1 J ll r ( (2 , + 1) 2_a) r ( ( 2 Z + 2)2-S)r(l-(2i + l)J + a)r(l-(2Z + 2)2 + a) 

(2.5.5) 

It is interesting to point out that the c 2 , C3, C4 terms in ?7 a , D a are related to the bound- 



ary if-matrix of the 6-vertex model[21| via the vertex-SOS intertwiner[p2], p0|, 23]. In 
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this light, one can regard the term c\ to be the solution special to the SOS/RSOS 
representation of the BYBE. In fact, it is not clear whether there exist an intertwiner 
that relate it to the K-matrix of the 6-vertex model. It is then natural to wonder if 
similar special solution occur in the elliptic case[p^]. 

Also intriguing is that restricting SOS to RSOS corresponds to 

C2 = C3 = C4 = 0, Co = ci sin(a7) . (2.5.6) 

3 Commuting transfer matrix 

Following the technique proposed in [Q] for the vertex model, one can similarly construct 
a family of commuting transfer matrix for the RSOS/SOS model with boundary. 

To start, it can be shown that if R% is a solution to the BYBE in the RSOS/SOS 
form then the following combination 

J2 St(u - Ul )S a c Ku + Ul )R a bc (u) (3.1) 

a 

also satisfies the BYBE, where (u) is the solution of the bulk YBE given in eqn. ( |2.1.2| ) 
and «i is an arbitrary parameter. The proof is essentially the same as the vertex case 
given in Q and we shall not repeat it here. 

It is convenient to think of the BYBE as the defining relation of some associative 
algebra generated by the symbol R bc . So the solutions given in the previous section 
correspond to particular representations of this algebra where the "quantum space" 
is trivial and the auxiliary space is the space of a one-step path V\ on a truncated 
Bratteli diagram with ab and ac being in- and out-state, respectively. In this context, 
the above "decorated" solution then corresponds to a representation whose quantum 
space is isomorphic to V\ that is formed by the nodes /, b (d, c). Clearly, the above 
construction can be repeated for an arbitrary number of times (say N + 1) giving a 
boundary i?-matrix that acts on Vn, the collection of iV-step paths on a truncated 
Bratteli diagram. We shall denote such a solution as R^ c which should be regarded as 
an operator acting on Vn- Its matrix element is explicitly given by 

N , 

Rt(u) « « = 5 , 5 baN 5 <WY [S a ' a * , (u-m) 

x SVI^ + OW" "(«) ( 3 - 2 ) 

a i-l a i J a o a o 

which has the graphical representation given in Fig.(Q). It carries N parameters U{ 
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Figure 4: "Decorated" boundary scattering matrix 
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from the bulk S-matrices and additional ones from R^ c . To form a commuting transfer 
matrix out of R^ c , like in the vertex case, one has to combine it with another BYBE 
solution (denoted here as R) as follows 

T(u) „ » = Y R%(u)Ri Ju) „ „. (3.3) 

a,b,c 

Hence the transfer matrix T(u) is again an operator acting on Vn- 
To show that 

[T(u),T(u')]=0, 
one inserts four bulk 5-matrices using the unitarity condition 

Y,S%(u +u)S%(-u -u)^8 ac , 

a 

and the crossing-unitarity condition 

S%(u + u)S%(r, - u + 7 - u) oc S M , 

a 

into T(n)T(u ). Then, one uses the BYBE to permute the R's, and the -R's. Because 
of the argument 7 — u + 7 — u that appears in the crossing-unitarity condition, one can 
take 

RUn) = RUl-u). (3.4) 
4 Diagonalization of the Transfer matrix 

So far, we managed to obtain solutions to the BYBE and construct the corresponding 
commuting transfer matrix. It would be necessary to diagonalize the transfer matrix 
in order to study the statistical models given by these solutions. For applications to 
field theory, diagonalization of the transfer matrix is also needed in order to write down 
the thermodynamic Bethe ansatz equation. For this purpose, a systematic approach 
generalizing the algebraic Bethe ansatz for the case of the periodic boundary condition 
has been devised in However, the method relies upon the conservation of the S z 
in the vertex language and is thus applicable only to diagonal boundary scattering 
theories. Therefore, we shall consider only the diagonal scattering solution and adapt 
the algebraic Bethe ansatz method devised in ||, along the line of |jlS|l , to the SOS 
model with boundary. The algebraic Bethe ansatz relies upon the existence of some 
pseudo- vacuum, which in the vertex model, is a state with either all spins equal to 1 /2 
or —1/2. In the SOS model, such a state corresponds to a path which takes the form 
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of a 45°-oriented straight line in the Bratteli diagram. It is obvious that some heights 
on such a path, for lattice size N large enough (> p), have to exceed the bounds 1 
or p — 1. Therefore, this pseudo- vacuum does not belong to the truncated Bratteli 
diagram and the algebraic Bethe ansatz that we are going to use is applicable to the 
SOS model only. For consistency, we have to assume that "f/n is an irrational number. 



To diagonalize the transfer matrix given in eqn^^), we have to first write down the 



algebraic relations satisfied by R^ c . As before, we express the operator R£ c as follows 

(«j-«i-i)/ 2 7 



where 



Ri<«) = *»M|$ftJ 11 [$) VwJSM 

+8 bc (Sb,a+iUa(u) + (5 M _iD fl (n))) (4.1) 

N 

1Zn{u) = lZ{u) T\U{u — Ui)U(u + Ui) , 
i=i 

which, along with g-number factors, ensures the boundary crossing and unitarity sym- 
metry of R£ c . Here X£ c ,U a ,D a are now non-commutative operators that satisfy the 
algebraic relations encoded by the BYBE. (see appendix) Because R^ c is diagonal, the 
commuting transfer matrix can be written as 

N-l /r "i\ ( u J-«i+i)/ 2 7 / r._ 11 
T(u) = K(j- U )K N (u) J] [Li^lT-ujUwH 

+l±p-D b+1 ( 1 -u)r> b+1 (u) S j . (4.2) 

It is important to bear in mind that the RSOS heights ao, ag on the "bottom" boundary 
(see Fig.(||)) are set to be the same since we consider only diagonal scattering. Let us 
denote them as a, while the heights on the other boundary are taken to be b as evident 
from the above equation. Hence the statistical model associated with the transfer 
matrix is the one defined on a square lattice whose boundary heights are shown in 
Fig.(||), where heights denoted by open-circles are summed over. 

Consider the column of heights in Fig. (|j) where aj = a + j and denote this state 
as u>a +N ■ By construction, uj^ +n is an eigenstate of XJ a+ N-i{u) and D a+ 7v+i(^) since 
the top and bottom heights of U a+ 7v_i(u)u;° +Ar and Y) a+ N + i(u)ui^ +N are a + N and 
a respectively. Furthermore, it is annihilated by a+N^ u ) ^ ue ^° ^ ne constraint 

that neighboring heights differ by ±1. Explicitly, 

U a+7V _ 1 (uK+ JV = UX(u)u*+ N 
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Figure 5: Lattice generated by the transfer matrix 



Y a+JV+l 
^a+N+2,a+N^ 



B a+N+1 (u)io a a +N 

where the eigenvalues are given by 
N ( N 



with 



v N Au)< +N 





0-1 



A a +i{ui) 
B a +i{ui) 



II Aa+iiUi) B a+j (Uj) Yl C a+k{Uk) \ D a+ i(u) 
j=l \i=j+l J \k=l 

N 

+ Yl Aa+j(uj)U a -l(u) 
3=1 

N 

Yl C a+j (ui)D a+1 (u) 



sin((a + sin((q + i — 2)7) sin(tt — m) sin(n + m) 

sin 2 ((a + i — 1)7) 
sin 2 7 sin((a + i — 1)7 + u — Uj) sin((a + i — 1)7 + u + Uj) 
sin 2 ((a + i — 1)7) 
C a +i(ui) = sin(7 - u + tij) sin(7 - u - Ui) 



and U a -i(u),D a+1 (u) are given by eqns.( |2.4.19| ),( |2.4.20D . 
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The state X"+j^I 2 a +N(^) UJ a +N 1S non-zero, which corresponds to a column of spins 
where the "bottom" and "top" spins have heights a and a + N — 2 respectively. So to 
obtain a state whose top spin has height given by b, one has to act on lo® +N successively 
with M = (N + a - b)/2 operators X^I*_ l a+7V __ - +1 (Aj) , % = 1, ■ ■ ■ , M. We shall 
denote such a state as 



<(A) = X^. 2 (A 



a+N-1 



a+N 



(4.5) 



^a+N-2,a+N\ 

which is the Bethe ansatz state, and the set of parameters A = (Ai, ■ ■ ■ , Am) have 
to satisfy some consistency condition necessary for ip b a {\) to be an eigenstate of the 
transfer matrix. Notice that M is always an integer fixed by the heights a, b and N. 

The algebraic relations among X° ±l tfFl ,U a ,D necessary for our purpose are (see 
Appendix) 

Va-^X^u) = g la (n\u)X:^(,)U a+1 (;)+ 52a (;,n)X^ 2 ( U )U fl+1 (,) 
+ ff3a (u', U )X^ 2 (n' )D a+3 H + g4 a {u ,u)X a + a 1 +2 {u)B a+3 (u) 

D a+1 (u')X^ 2 («) = f la (u',u)X. a a + 1 +2 (u)B a+1 (u) + f 2a (u 1 u)X. a a + 1 +2 (u')n a+1 (u) 
+ / 3a (u',«)X^ 2 (n')U 0+1 H, (4.6) 



where 



gia(u ,U 
92a{u ,U 
93a(u',U 
94a(u ,U 
fla(u',U 
ha{u ,U 
ha (u',U 



sin(a7) sin((a + 1)7) sin(7 — u + u) sin(27 



u 



sin((a — 1)7) sin((a + 2)7) sin(7 — u — u) sin(u' — u) 
sin 7 sin((a + 1)7) sin(a7 + u — u) sin(27 — u — u) 
sin((a — 1)7) sin((a + 2)7) sin(7 — u — u) sin(n' — u) 
sin7sin((a + 1)7) sin(a7 + u + u) sin(27 — u + u) 
sin((a — 1)7) sin((a + 2)7) sin(7 — u — u) sin(u' — u) 
sin 7 sin(27) sin(a7 + u + u) sin((a + 1)7 + u — u) 
sin((a — 1)7) sin((a + 2)7) sin(7 — u — u) sin(n' — u) 
sin(n' + u) sin(7 + u — u) 
sin(7 — u — u) sm(u — u) 
sin7sin(n' + u) sin((a + 2)7 — u + u) 
sin((a + 2)7) sin(7 — u — u) sin(u' — u) 
sin7sin((a + 2)7 — u — u) 
sin((a + 2)7) sin(7 — u — u) 



The action of the transfer matrix on the state ^q(A) can be evaluated by commuting 
Ub_i(«) and D 6+1 (u) through X»+^I*_ 

-i j0+ jv+i-i(^i)- However, the presence of the 
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34a (u i u) term in the first commutation relation complicates matters considerably. Like 
the vertex case, it is desirable to define a new operator 



with 



Uo-i(m) = U _i(u) + h a (u)D a+1 (u) 
, sin7 sin((a — 1)7 + 2u) 



(4.7) 



sin((a — 1)7) sin(7 — 2u) 
so that the operators U a -i(u), D a+ i(u) satisfy simplified relations 

U a _i(ti')X2jj. 2 («) = aia(u\u)X%+\ 2 (u)tla+i(u) + a 2a (u ,u)X^ a 1 +2 (u')ij a+1 (u) 

+ a 3a (u\u)Xl+ a 1 +2 (u')D a+3 (u) (4.8) 
(3 la (u , u)X£& 2 (u)D a+1 («' ) + (3 2a (u , n)X^ 2 («' )D a+1 (u) 

(4.9) 



D a+1 (u')X^ 2 («) 



where 



a la (u ,u 
a 2a (u',u 
a 3a (u , it 

Pi a (u',U 
@2a{u ,U 
P3a(u ,U 



+ / g 3a («',«)x^ 2 (u')u 0+ i(«) , 

sin(a7) sin((a + 1)7) sin(7 — it' + it) sin(27 — u — u) 
sin((a — 1)7) sin ((a + 2)7) sin(7 — u — u) sin(u' — it) 
sin7sin(a7) sin((a + 1)7 + u — it) sin(27 — 2u) 

sin((a — 1)7) sin((a + 2)7) sin(7 — 2u) sin(it' — it) 
sin7sin(a7) sin(2n) sin(a7 + u + it) sin(27 — 2u ) 
sin((a — 1)7) sin((a + 1)7) sin(7 — u — u) sin(7 — 2u) 
sin(u' + it) sin(7 + u — it) 
sin(7 — u — u) sin(u' — u) 

sin7sin(2ti) sin((a + 1)7 — u + it) 



sin((a + 1)7) sin(7 — 2u) sin(it' — it) 
sin 7 sin((a + 2)7 — u — it) 
sin((a + 2)7) sin(7 — u — u) 



Moreover the eigenvalue (u) of U a _i(u) on lo^ +N ls much simplified and given 



by 

Z4» 



N 



sin(a7) sin(2it) sin(£ a + (a - 1)7 + it) sin(£ a + 7 — u) 
sin ((a - 1)7) sin (7 - 2u) sin(£ a + 07 + 11) sin(£ a - u) 



(4.10) 
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where use has been made of the relation 

B a+i {ui) + h a+i (u)C a+i (ui) = h a+ i-\{u)Aa+i{ui). (4.11) 
Another important relation we need is 

X«+ 1 +2 («)X«+| a+4 (u) = X£+V 2 («')X«+| a+4 («) (4.12) 

which implies that the Bethe ansatz state ij)\{X) is a symmetric function in the Aj's, 
which is useful in obtaining the Bethe ansatz equation. 

It is now straight forward to compute the action of U&_x(it), D& + i(ti) on ^„(A) 
giving 



U 6 _i(«)^(A) 
+ 



= ai6(it, Ai)ai6+2(ti, A 2 ) • • • a la+N - 2 (u, X M )U^ (u)ip^(X) 

(X2b(u, Al)ai6+2(Al, A 2 ) • • ■ QH +JV-2(Al, X\i)Ua (Al) 

+ a3b(n,Ai)^ib+2(Ai,A 2 )---^ia+Ar-2(Ai,AM)^(Ai) ^(A;Ai) 



+ 
+ 



«2b("u> Am )ai6+2(AM) A 2 ) • • • aia+Jv-2(AM, Am-i)W q (Am) 



+ «3&(^; Am)/3i6+2(Am, A 2 ) • • • Pio,+n-2(Xm, Am-iR (Am) ^(AjAm) 

(4.13) 



and 



D 6+1 (u)^(A) = Ai)/?i fe+2 (n, A 2 ) • • • Pi a +N-2(u, X M )V^ {u)ip h a (X) 

+ P?b(u, Ai)/3i fe+2 (Ai, A 2 ) • • • /3i a +iv-2(Ai, X M )T>% (Ai) 

+ /5 3 ; ) ('u,Ai)ai fe+ 2(Ai,A 2 )---ai a+ Ar_2(Ai,AAf)^i v (Ai) Va(A;Ai) 
+ ••• 

+ #26 («j Am)#i6+2(Am) A 2 ) • • • /?i a+ Ar_ 2 (AM, Am-i)£> (Am) 

+ /33fe( - ",AM)ai6+2(AM,A 2 )---ai a+ Ar_2(AM,AM-l)^i V (AM) V'o(^)^iw) > 

(4.14) 

where the state ^(A;Aj) is defined as in eqn.( |4.5p with Aj replaced by u. Combining 
the above in the transfer matrix (see eqn.(f4.2|)), 



(4.15) 



sin(27 - 2u) sin(£ fc - u) sin(£ fc + 67 + ii)A>+i (7 - u) , . 



+ 



sin(7 — 2u) sin(^ + 7 — n) sin(£{, + (6 — 1)7 + u) 
sin((6- l)7)J7 6 _i(7 - u) ~ 



sin(&7) 



-U 6 _i(«) 
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the state ^(A) is an eigenstate of the transfer matrix if the coefficients of ^„(A; \)i * = 
1, • • • , M vanish, which gives the following Bethe ansatz equation to be satisfied for Aj's 

sin(g fe - 7 + Xj) sin(g b + (b + 1)7 - Xj) sin(g a + 7 - Xj) sin(£ a + (a - 1)7 + Aj) 
sin(£ b - Aj) sin(£ b + 67 + Aj) sin(£ a + Aj) sin(£ a + 07 - Aj) 

sin(Aj + Aj) sin(7 + Xj - Xj) sin(7 - Xj + u k ) sin(7 - A— life) 

J-- 1 - sin(27 - \ - Aj) sin(7 - Ai + Aj) M sin(Aj - u fc ) sin(Aj + u k ) 
j = 1 

for i = 1, •••,M 

and the eigenvalue A* (u, u; A) of the transfer matrix is given by 

A b a (u,u;X) = 7^(7 - u)1l N (u)U b ^i{~{ - u)D a+l (u) 

sin(27 — 2u) sin(^ — u) sin(^5 + 67 + 14) 



(4.16) 



,sin(7 - 2u) sin(£{, — 7 + u) sin(£ fe + (b + 1)7 - u) 

n A sin(u + Xi) sin(7 + u - Xi) SL 
^ gin(At _, )sin(7 _,_ At) H ^(7 -« + «*) -n(7 - « - «*) 

sin(2u) sin(£ a + 7 - n) sin(£ a + (a - 1)7 + u) 
sin(7 - 2u) sin(£ a + u) sin(£ a + 07 - u) 

^ sin(u + A« - 2 7 ) sin(7 - u + Xi) . \ 

x I I —, \ \ • / — 1 1 sm ( u ~ u k) sm ( n + Uk) ) 

iJ{ sm(u- Xi) sm(j - u - Xi) fc A =1 J 

(4.17) 

In the above analysis, £ a and ^ are respectively the free parameters associated with 
the boundary i?-matrices at the bottom and top boundaries. They need not be related 
at all, hence more generally one should write them as £<T and to distinguish their 
origins. 

5 Discussion and open problems 

In this paper, we present the general trigonometric RSOS/SOS solution to the BYBE. 
By comparing them with the corresponding solution in the vertex representation, one 
may be able to obtain useful information about the vertex-SOS transformation matrix 
for the BYBE ^0|. Indeed, for the diagonal solution, in the limit £ a — > ±ioo, U a -\ 
and D a+ \ are equal so they contribute as an overall factor for the transfer matrix, 
since each bulk weight is invariant under the action of the quantum group XJ g sl(2) 
symmetry (q = e* 7 ), the transfer matrix possesses the quantum symmetry. In this limit, 



the vertex-SOS transformation [25] is well known and has a precise meaning in terms 
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of the Clebsch-Gordan coefficients of U 9 sZ(2). Thus for generic £ a , the vertex-SOS 
transformation can be considered as an extension of that of the \J q sl(2) case. 

The SOS/RSOS solutions indicate that the diagonal solution is not a special case 
of non-diagonal solution, in particular, there is no way of adjusting the free parameters 
to make all the off-diagonal scattering weights Xff vanish. The diagonal solution is the 
most favorable case to be studied as we have demonstrated how the transfer matrix 
of the SOS lattice built up from this solution may be diagonalized. On the other 
hand, it is not obvious how such method can be applied to the RSOS case. In the 
limit 7r/7 becomes a rational number, solutions to the Bethe ansatz equation for the 
SOS model should contain those for the RSOS. However, except at the special point 
where there is a V q sl(2) symmetry [26|, it is not clear how the RSOS solutions may 
be extracted. It would be an interesting challenge to extend the idea to any £ a . The 
RSOS model is indeed a very interesting case to consider; it has been shown that for a 
different geometry where bulk faces are oriented at an angle 45° with respective to the 
boundaries, the RSOS boundary condition where all spins have height a corresponds in 
the continuum limit to the boundary conformal state \h\ A ) and that with all boundary 
spins and their neighbors have respectively heights a,a + 1 corresponds to \h a \) 



The lattice model that we considered has a different gemeotry from that in [27], however, 
we believe that the difference is not significant in the scaling limit. Then the former 
boundary condition in fact can be obtained as the £ a — > ±ioo limit of the RSOS 
diagonal solution since the weights U a -\, D a+ \ become the same. While the latter 
boundary condition can be obtained with £ a = u. So for generic £ a , the RSOS diagonal 
solution is in fact a mixture of the two above-mentioned boundary conditions and it 
would be interesting to examine which boundary conformal state it corresponds to in 
the scaling limit. Similarly, the non-diagonal solution, which gives to some extent a 
free boundary- like condition, may correspond to boundary conformal state \h r>s ) with 
r, s ^ 1 |2£| in the scaling limit. 

As a scattering theory, our results should be related to the conformal boundary 
conditions of the boundary minimal CFT and perturbations by relevant operators in 
the bulk and on the boundary. The non-diagonal solution, which has no free param- 
eter, can be interpreted as the bulk $13 perturbation of the CFT with free boundary 
condition where all possible spins are allowed on the boundary. Additional integrable 
boundary perturbations introduce CDD-like factor in the scattering amplitudes. Each 
diagonal solution describes a perturbed CFT with a fixed boundary condition where the 
boundary B a has fixed spin a. The parameter £ a in the solution should be related to the 
coupling constant of the boundary perturbing field. It is an open problem to relate our 
solutions to the specific conformal boundary conditions and boundary perturbations of 
generic minimal CFTs. The case p = 4 has been analyzed in f29|| . 

One can also generalize our results to the coset CFTs SU(2)fc (g) SU(2)//SU(2)fc + j 
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perturbed by the least relevant operator in the bulk and by some boundary fields. The 
bulk-scattering matrices are given by 

S = <SRSOS(fc+2) ® <Sr,SOS(Z+2) , 

where SrsosO) is the RSOS S-matrix of the kinks. In this theory, particles carry two 
sets of RSOS spins and can be represented as \K a b) \K a > b r) [12|. For the BYBE 
eqn.( [2.1."8| ) with the above bulk scattering matrix, the boundary i?-matrix given by 

R = ^RSOS(fc+2) ® -RrSOS(Z+2) ; 

is a solution, where -Rrsos(p) is the i?-matrix given in eqn.( [2.1."9l) . In particular, with 
k = 2, this is the N = 1 super CFTs with boundary perturbed by the least relevant 
operator and the i?-matrix is given by that of the tri-critical Ising model tensored with 
that of the RSOS. 

To answer some of the questions raised, finding the Bethe ansatz equation and 
diagonalizing the transfer matrix for the RSOS case are the essential first step. It is 
quite likely that one needs an alternative method such as functional approach for this 
purpose. 
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Note added in proof 

After completing this work, we learned that some of the results presented in section 2, 
3 have also been independently obtained in |hJ ^4|] as their trigonometric limit. 



Appendix 

We present here the complete set of algebraic relations satisfied by the operators X^ c («), 
TJ a (u), D a (u). These relations are obtained from the BYBE using eqn. (|4.l| ) for the 
boundary scattering matrix. Expanding the BYBE equation and considering the vari- 
ous allowed heights, we get 

ci«(u + )ci (u_)U a _iUa_ 1 + ci a (w_)X^ X^i 2 + c 3a (u',ii)U a _iD' a+1 
= ci a (u + )ci a (n_)U| 1 _ 1 U a _i + ci a (u_)X a a _- 2 ^X^-?_ 2 + c 3a (u , u)D^ +1 U a _i 
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C2a(u + )c2a(U-)~Da+l~D'a + l + c 2a )X^ a X^2 + C 3a (u , MjD^iU^j 
= C 2 a(n+)c2a(n_)D^ +1 D a+ i + C2a(u_)Xj+_ 1 L 2 X at2,a + C3a(«' , u)\j' a _{D a+1 

C la /_( U+ )U a _ 1 Ul_ 1 + C 2a / + (n_)U a _ 1 Dl +1 + f-K7a-2^a-2,a 

= c 2a /_(n + )Dl +1 D a+1 + c la / + ( U _)Ul_ 1 D a+1 + f-X'^Ktla 

c 2a f-(u + )T> a+1 D' a+1 + Cla / + (n_)D a+1 Ul_ 1 + /-X:+| 2 XX 
= c la /_(n + )Ul_ 1 U a _ 1 + c 2a / + ( U _)Dl +1 U a+1 + f_X' a a -} 2 X a a Zi a 

c 2a _ 2 (n + )X^ a 1 _ 2 X^ a + Va-iK-i = c 2a _ 2 ( u+ )xJ a L 1 2 x^: 2>a 

+u a _ 1 u;_ 1 

cia+2(^+)X^ 2 X a a + 2 : | a + D a+ iD a+1 = ci a+2 (n + )X a ^ 2 X^+ 2!a 

+ Dl+lDa+l 

va -v"'a+2 "V" /a -y-a+2 

■^o-l,o+l^o+l,o+3 — -^0-1,0+1^0+1,0+3 

X^^a-l + C la («+)U a _iX^ 2 = Cia^-JX^Da-x 

+ ci a (u+)ci a (n_)U a _ 1 X^ a _ 2 + c 3a (u , «) 

X^ a+2 Ul +1 + c 2a (u + )D 0+1 X^ 2 = c 2a (u_)X^ 2 U a+1 
+ c 2a (n + )c 2a (n_)D; +1 X^ 2 + c Za (u ,u)\j' a _{X a + a \ 2 

C2a(u-)Uo+iX£g a + C2a(n+)c 2a (n^)X^ )a Dl +1 
+ c 3o (</, u)XlXi a V' a -i = K+lKXla + c 2a (u+)X^ D a+1 

c la (u_)D a _iX^ a L 1 2 + ci a (n+)ci a (n_)X^ 1 _ 2 Ul_ 1 

+ ^(u'^x^X+i = d^x-:^ + Cl0 (« + )x;°_- 2 ; a u _i 

/+Ua-iX^ 2 = /_X^ 2 U a+1 + c la (u^f + V' a _ 1 Xl+ a 1 +2 
+ c 2o (u + )/_D; +1 X^ 2 

Z+Da+iX^a = /_X^ 2 D a _! + c^^D^X^ 1 ., 

+c la ( U+ )/_u:_ 1 x^ i _ 2 

f-U a+1 X' a %) a + ^(u-tf+X^X-i + c 2a (u + )f-XlXl a T>' a+1 

= /+ X a+2la U a-l 
= /+ X a-2,a D a+l 
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where 



c\ a {u) 

C2a{u) 

c 3a (u',u) 

f± 

u± 

Here again, we abbreviate 



sin7sin(a7 — it) 
sin(a7) sin(7 — it) 
sin7sin(a7 + it) 
sin(a7) sin (7 — u) 

sin((a — 1)7) sin((a + 1)7) sin(n + u) sin(it — it) 



sin 2 (07) sin(7 
sin(n ± u) 
sin(7 — it' =F u) 
u ±u . 



it) sin(7 — it + u) 



u: 



U a («) 
U a (u) 



and similarly for D a , X^ c . 

Among them eqns. (l5.9| ),( |5,15 ) are of interest to us, with the help of the latter the 



former can be turned into the first equation in eqn.(4.£) which is more convenient for 
the algebraic Bethe ansatz computation. 
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